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Abstract. Solutions ol the classical 4 -theory in Minkowski space-time are ana- 
lyzed in a perturbation expansion in the nonlinearity. Using the language of Feynman 
diagrams, the solution of the Cauchy problem is expressed in terms of tree diagrams 
which involve the retarded Green's function and have one outgoing leg. In order to 
obtain general tree diagrams, we set up a "classical measurement process" in which 
a virtual observer of a scattering experiment modifies the field and detects suitable 
energy differences. By adding a classical stochastic background field, we even ob- 
tain all loop diagrams. The expansions are compared with the standard Feynman 
diagrams of the corresponding quantum field theory. 
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1. Introduction 

The dynamics of interacting quantum fields is most successfully described pertur- 
batively in terms of Feynman diagrams. One distinguishes between tree diagrams, 
which are finite, and loop diagrams, which typically diverge. Performing a suitable 
"classical limit" h — > 0, the loop diagrams drop out, and only the tree diagrams re- 
main. Therefore, it is a common opinion that tree diagrams are inherent already in 
the classical theory, whereas loop diagrams describe the particular effects of quan- 
tized fields. However, no hints can be found in the literature on how this statement 
could be made precise. This was our motivation for approaching the question coming 
from the classical theory: We consider a classical field, being a solution of a nonlinear 
hyperbolic partial differential equation. Treating the nonlinearity perturbatively, we 
obtain an expansion which we express in the language of Feynman diagrams. Then 
the question is, does the resulting expansion contain all tree diagrams of quantum field 
theory? Are there differences in the tree expansions? Do loop diagrams appear? Is 
it really impossible to obtain loop diagrams within the realm of classical field theory? 
Our hope was that analyzing these questions would give a better understanding of 
what the essence and physical significance of field quantization is. In particular, we 
wanted to clarify how the nonlinear dynamics of classical fields and the corresponding 
nonlinear scattering theory fit into the framework of a linear dynamics of quantum 
fields on the Fock space. 

Our results are unexpected in several ways. First, rewriting the solution of the 
Cauchy problem in terms of Feynman diagrams, we only obtain diagrams with one 
outgoing leg, in which all lines correspond to retarded Green's functions. Several out- 
going legs are obtained only by considering a scattering experiment and setting up a 
"classical measurement process" , in which the wave is modified with an inhomogeneity 
at large times, and the difference of the energies before and after this modification are 
measured. Furthermore, in order not to distinguish a direction of time, we replace 
the perturbation expansion involving the retarded Green's function by the so-called 
global perturbation expansion. This gives complete agreement with the tree diagrams 
of quantum field theory, except that, instead of the Feynman propagator, other, nec- 
essarily real-valued, Green's functions appear. In this sense, we can confirm that tree 
diagrams are in fact inherent already in the classical theory. However, this statement 
is true only in the special setting of the classical measurement process, and only up to 
the differences in the choice of the propagators. 

As another surprising result, we find that it is indeed possible to obtain loop di- 
agrams in a purely classical context. To this end, we consider the classical field in 
a stochastic background, as described by a Gaussian ensemble of classical solutions. 
The resulting expansion in terms of Feynman diagrams includes all loop diagrams 
of quantum field theory. Our expansion differs from that of QFT in that the Feyn- 
man propagators are to be replaced by real-valued Green's functions and fundamental 
solutions, with a specific combinatorics. Moreover, there are differences in the combi- 
natorial factors for higher order loop diagrams, as we explain in a simple example. 

Working with a stochastic field has some similarity with the approaches to explain 
quantum effects by adding a stochastic term to the classical equations (see for example 
Nelson's stochastic mechanics [14] or [UHl]). However, in contrast to these approaches, 
we do not modify the classical equations. We merely add a stochastic background field 
to the initial data. This additional field participates in the classical interaction as 
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described by the classical field equations. The nonlinearity of the mutual interaction 
gives rise to the loop diagrams. The physical picture is that the macroscopic field 
is superimposed by microscopic fluctuations, which can be thought of as a classical 
analog of vacuum fluctuations. These microscopic fluctuations can be observed only 
indirectly via their influence on the dynamics of the macroscopic field. 

We point out that our model of a classical field theory in a stochastic background 
field is clearly not equivalent to a quantized field. In particular, our wave modes are 
classical, whereas in quantum field theory they correspond to quantum mechanical 
oscillators. But our point is that on the level of Feynman diagrams, these differences 
might not be visible. In particular, the radiative corrections might not necessarily be 
a test for the quantum nature of physical fields. 

As is common in perturbative QFT, our treatment is formal in the sense that we 
do not care about the divergences of diagrams and disregard all convergence issues 
of the expansions. Thus all our expansions are to be understood as formal power 
series involving symbolic expressions. Moreover, in order to keep the setting as simple 
as possible, we always restrict attention to the real massless ^-theory in (3 + 1)- 
dimensional Minkowski space. But our methods and results immediately extend to 
any other bosonic field theory in arbitrary dimension. 

The paper is organized as follows. In Section [2] we provide the necessary background 
on the Cauchy problem in classical field theory and on Feynman diagrams in QFT. In 
Section [3] we consider expansions of the classical field in terms of Feynman diagrams. 
Taking free fields as the starting point (Section I3.1f) . we express the solution of the 
Cauchy problem in terms of Feynman diagrams (Section I3.2() . We then set up the 
"classical measurement process" and derive an expansion in terms of tree diagrams 
(Section EMS]). 

Section [4] is devoted to the classical field theory in a stochastic background. After 
introducing the free stochastic background field (Section 14. ip . this field is included in 
the perturbation expansion to obtain loop diagrams (Section 14. 2ft . The comparison to 
the loop diagrams in QFT is given in Section [4.31 

Finally, in Section[5]we give an outlook on more realistic theories including fermions. 
We also outline potential applications of our methods to constructive field theory. 

2. Preliminaries on 4 -Theory 

We introduce classical (/> 4 -theory in the Lagrangian formulation. The Lagrangian L 
is given by 



where is a real- valued scalar field. Integrating the Lagrangian over Minkowski space- 
time gives the action S, 



We are working in units where c = 1 = H and denote a chosen length scale by I. Then 
the field has dimension of inverse length, whereas the coupling constant A and the 
classical action are dimensionless, 



S 




M = r 1 



and [A] = It 



\S\. 
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Considering critical points of the action, one obtains the Euler-Lagrange (EL) equa- 
tions 

□<A = "^ 3 (2-1) 
b 

(where □ = e)f — A]g3 is the wave operator). 

According to Noether's theorem, the symmetries of the Lagrangian correspond to 
conserved quantities. In particular, the symmetry under time translations gives rise 
to the conserved classical energy E, 

m = J (l^ + l i v ^i 2 + 1 ^) . ( 2 - 2 ) 

having the units 

[e] = r 1 . 

2.1. Basics on the Classical Cauchy Problem and Scattering Theory. In the 

Cauchy problem, one seeks for solutions of the nonlinear hyperbolic partial differential 
equation (|2.ip for initial data ((f), <f)\t given at some time to- It is a classical result 
that the Cauchy problem is well-posed, meaning that there is a unique solution for 
short times in a suitable Sobolev space. A general method for the proof is provided by 
the theory of symmetric hyperbolic systems (see |10t 121]), The theory of symmetric 
hyperbolic systems also reveals that the solutions of hyperbolic partial differential 
equations propagate with finite speed, showing that the physical principle of causality 
is respected. In general, the solution of the Cauchy problem will not exist for all 
times, as singularities may form. It has been a major topic of mathematical research 
to obtain global existence results and estimates of the solution for long times (see for 
example the textbooks [HI HE1 Q2]). Here we shall disregard all analytical questions by 
simply assuming that our solution exists for all times. 

In scattering theory, one is interested in the asymptotic behavior of the field in the 
limits t —> ±oo. Scattering theory has been developed mainly for linear equations 
(for good references on linear scattering theory in the mathematical physics literature 
see [TU [El [TTj, [2, 24J). In the standard setting, the dynamics is described by a linear 
evolution equation of the form of the Schrodinger equation 

id t iP = H(t)i>, (2.3) 

where H{t) is a self-adjoint linear operator on a Hilbert space (., .)). Then the 

corresponding time-evolution operator U(tx,t>i) defined by 

U(h,t 2 ) : %{t 2 ) -+ 5C(ti) : ^ V>k 

is a unitary operator between Hilbert spaces. One wants to compare the interacting 
dynamics (12. 3D with a free dynamics of the form 

id t ^ = H (t)^, (2.4) 

where Hq is a self-adjoint linear operator on a Hilbert space (^Kq, (., .)). Usually, the 
operator Hq is so simple that the time evolution operator f7o(ti, t 2 ) is explicitly known. 
Then the wave operators are introduced by considering the late-time limits 

OtV>:= lim U(t ,t) U (t,t )i/> : JCo -> Oi = JC(t ) , (2.5) 

t— >±oo 

where for large times we identified the Hilbert spaces !Ko and %(t). Then the scattering 
operator S is defined by 

5 = (0_)*^+ : IKq -> IKq ; (2.6) 
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it maps the free outgoing field to the corresponding free incoming field. 

As the -theory is non-linear, the formalism of linear scattering does not imme- 
diately apply. In order get a rigorous connection, one would need to prove that the 
wave dissipates for large times, going over to a solution of the free wave equation. This 
requires decay estimates in the spirit of [HI [6] ; see also [T71 Chapter 9] and [201 Sec- 
tion 6]. In this paper, we shall not enter dissipation and decay estimates, but instead 
we will restrict attention to a systematic treatment of nonlinear classical scattering on 
the level of formal perturbation expansions. We will enter this analysis in Section [3l 
In order to get into the position to compare the different settings, we now briefly 
summarize the basics on quantum scattering. 

2.2. Basics on Quantum Field Theory and Feynman Diagrams. We now re- 
view the basics on the corresponding (^-quantum field theory. For brevity, we keep 
the presentation on the level of formal calculations (for attempts for making these 
calculations mathematically sound see for example [7])- We begin with the formalism 
of canonical quantization. To quantize the free field, we first write the general solution 
of the classical equation D0 = as a Fourier integral, 

(2 - 7) 

where we set oj = \k\. The fact that our field is real means that the Fourier coefficients 
satisfy the relation 

M",k) = M-",-k)- (2.8) 
For the quantization, the Fourier coefficients are replaced by operator-valued distribu- 
tions according to 

4>o(uj, k) — > 2ir V%uj a(k) , <Po{— —k) — > 2ir V2uj a){k) , 

where a(k) and a'(k) satisfy the commutation relations 

[a(k),a\q)] = (2<k) 3 5 3 (k - q) , [a(k),a(q)] = = [a\k), a\q)} . (2.9) 

We denote the resulting quantized field by <f>o, 

k{x) = / tSs ^ ( a W + flt (^) ■ ( 2 - 10 ) 

As an immediate consequence of (|2.9|) and (|2.1(jp . the quantized field has the following 
properties: 

□0o(a;) = (2.11) 
[Mt, x), Q (i, y))=0= [dtMt, x),d t 4>o(t, y)] (2.12) 
[d t Mt,x),Mt,y)] = -i5 3 (x-y) . (2.13) 

Noting that dt<ft is the canonical momentum, equation (|2.13p is the analog of the com- 
mutation relation [p, q] = —i of quantum mechanics. For general space-time points x 
and y, a short calculation yields the commutation relations 

[Mx),My)}=^K (x-y), (2.14) 
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where Kq(x) is the distribution 

i r d 3 k 



— ILUt+ikx iuit — ikx \ /n Irt 

(2vr)4 J 2u \ 6 6 ) ' [ b) 

We remark that the equal-time commutation relations f)2. 12|) and (|2.13p immediately 
follow from ()2.14j) and (|2. 16|) by differentiating with respect to time and taking the 
limit x° — y° — > 0. The Fourier representation (|2.15p shows that Kq(x) is a distribu- 
tional solution of the wave equation. Moreover, this distribution is causal in the sense 
that it vanishes for spacelike x. This becomes apparent from the representation 

where S v and S A are the advanced and retarded Green's function defined by 

1 ( e —iut+ikx _ e iuit+ikx\ \[ I ~ (I 



(2vr)3 2oj V ) (2.18) 

if t < 

5 v (x) = S A (-x) . (2.19) 

Usually, the operator algebra which is generated by a(k) and a'(k) is realized on a 
Fock space (3~, (.|.)). This is a Hilbert space with a distinguished normalized vector |0) 
(the "vacuum state") having the following properties: 

(i) For all k, the operator a(k) annihilates the vacuum: a(k)\0) = 0. 

(ii) The vacuum is cyclic in the sense that by iteratively applying the opera- 
tors a'(k) to |0) one generates a dense subspace of 3". 

(iii) For all k, the operator a'(k) is the formal adjoint of a(k) with respect to the 
scalar product (.|.). 

The properties (i)-(iii) together with the commutation relations (]2.9p determine the 
scalar product. For example, setting 

\k) := a\k)\Q) and (k\ := (\k)\ . ) , 

we obtain 

(k\q) = <0|a(£)at(g)|0> = <0| [a(k), a\q)\ |o) 

= (2^) 3 5 3 (k - q) (0|0) = (2vr) 3 5 3 (k-q). 
Thus, introducing for a Schwartz function / the vector 

/J 3 k 1 -> -> 

one obtains 

1 f d 3 k 
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In order to introduce the interaction, one wants to construct new field operators eft 
on J which solve the nonlinear classical field equation, 

n$(x) = -\0(xf . (2.22) 

This can be accomplished by the ansatz 

<j>(t,x) = U(tf 0o (t,x)U(t) , (2.23) 

where U(t) is a formal solution of the equation 

id t U(t) = -H int (t) U(t) with H int (t) :=jf Mt,xf d z x . (2.24) 

J t=const 

Indeed, applying the wave operator to (|2.23[) and using (|2.24p together with the com- 
mutation relations (|2. 12j) and (|2. 13j) . a straightforward calculation shows that (|2.22|) 
is satisfied (for details see for example [22| Section 3.5] in the similar context of a 
pion- field or [151 Chapter 4]). 

Writing the interaction in the form ()2.23f) and ()2.24p is particularly useful for describ- 
ing a quantum scattering process. To this end, we assume that in the limits t — > ±oo, 
the dynamics goes over to that of free fields. Thus the asymptotic states are vectors 
of the free Fock space (3~, (.|.)). We denote the free field operators in the asymptotic 
regions by <j)\ n and cp out (note that they are solutions of the free field equation (|2.1ip ). 
According to common notation, the Fock space (3~, (.|.)) with the basis generated by a\ n 

and a£ ut is denoted by (3~i n , (.|.)i n ) and (S^utj (-l-)out)) respectively. By construction, 
these two Fock spaces are identical, but they come with different orthonormal bases. 
In order to describe the transformation of one basis to the other, one considers the 
solution of the above equations (12.231) and (12.241) with 4>q = <p- m and 

lim U(t) = 1 . (2.25) 

t— >— oo 

Then U(t) can be given in terms of a time-ordered exponential, 

U(t) = Texp (i f H- mt (r) dr) = Texp ( ^ j ^ ^ n (y) A d A y\ 




d A y n T (0 in (yi) 4 ••• ^m(yn) 4 ) 
^o n! ^ 4! ^ J {y°<t} J{i&<t} V ' 

The scattering operator S is defined by 



S := lim U(t) (2.26) 

t— >oo 

= E ^ (ji) J d ^ ■ ■ ■ J ^ T {^vi) A ■ ■ ■ ^n(yn) 4 ) • (2.27) 

Then, according to (|2.23|> . 

<Pout = in S . 

Moreover, noting that H[ nt is a symmetric operator on the Fock space, one sees that S 
is a unitary operator on 3~, 

S : 3" — >• 3~ unitary 
| out) i y | in) . 
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We remark that this construction is very similar to the procedure in classical linear 
scattering theory. While the scattering operator in the linear classical setting f)2.6|) is 
a linear operator on the Hilbert space "Kq of classical fields, the quantum scattering 
operator (|2.26p is a linear operator on the Fock space 3~. Since here we are working in 
the interaction picture, the free dynamics Uq is the identity. Then the relations (|2.5p 
and (|2.6p reduce to S = U(— 00,00). This is consistent with (|2,26p if in view of (|2.25p 
we identify U(t) with U(—oo,t). 

The scattering amplitude of a scattering process can be expressed as a matrix ele- 
ment of the scattering operator with respect to the in-basis, 

(A>ut I a in ) = (An |5| Ofa,) (2.28) 



^ n\ V 4! 

n=0 



d 4 yi ■■■ d 4 y n //3in T ^ in (yi) 4 ■ ■ ■ <fci(yn) 4 ) a in \ 



= (An I a in ) - ^ Jd 4 y (An I 0in(y) 4 | a in ) + 0(A 2 ). 

Generating the in-basis by iteratively applying (f)- m to the vacuum (LSZ reduction 
formalism; see for example [E]), all expectation values can be expressed in terms of 
the n-point functions defined by 

G (XI, . . . ,X n ) = (O|T(0 in (xi) • • • 4>in(x n ))\0) ■ 

Using the commutation relations, the n-point functions can be written as sums of 
products of the two-point functions (Wick's theorem). A short computation using (|2.9p 
and (|2.10p yields for the two-point function 

G {x,y) = iA F (x-y), (2.29) 

where is the Feynman propagator, 

d 4 k e~ ikx 



A_p(x) := lim 



(2.30) 



e\o J (2vr) 4 k 2 + is 

(2^3 2^ ( 0(t) e ~ iU)t+ikx + (-*) eWt ~ i] 

The Feynman propagator satisfies the defining equation of a Green's function 

-DA F (x) = 5 A (x) . 

However, in contrast to the advanced and retarded Green's function, it is complex- 
valued, and it is non-zero in space-like directions. 

For the systematic treatment of perturbation theory, it is most convenient to work 
with a generating functional in the path integral formalism. We thus introduce the 
generating functional 

Z x [j] = J e ^S[4>]+r4>) where j ■())■= j j(x) <j){x) d A x . 

The generating functional is related to the vacuum expectation value of canonical 
quantization by (see for example [23] and [8]) 

Z x [j] = (0 j Texp (~ 1 0f» d A x + i f j(x) 4>in(x) d 4 x) j 0) . 
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In order to compute Z\\j], we rewrite the generating functional as 

Z X [j] = exp Jd*x(- z^) 4 ) Z [j] , (2.31) 
where Zq[j\ is the generating functional of the free quantum field theory: 

Z \j] -=1^1 V<p exp (i J (~\<t>{x) ((a - ie)<P(x)) + j{x) 4>{x))(fx 

= lim / T>/j, e ((j)) e^m+r^) _ (2.32) 

e\0 J 

Here Sq is the free action, 

SM = \ J(d^)(x) (d^)(x)d 4 x, 
and T>u e (4>) is the Gaussian measure 

This path integral can be computed explicitly to obtain 

Z [j] = exp {- % - J j j(x) A F (x - y) j(y) d 4 xd 4 y^j (2.33) 

(we use the convention Zq[0] = 1). In particular, 

S S 

G (x,y) = - Z [j] =iA F (x-y), 

in agreement with (|2.29p . Wick's theorem can be expressed as 



G { Xu ...,x n ) = {-*T~-~m\^ 



6 6 e -Ufj(x)A F (x-y)j(y)d i xd i y\ 



(2.34) 



showing that the n-point functions can indeed be expressed as sums of products of the 
Feynman propagator. 

The scattering amplitudes are most conveniently expressed in terms of the interact- 
ing n-point functions G\, which are obtained by taking functional derivatives of the 
generating functional with interaction, 

Gx(xu ... iXn y.= ( _ i) n^..._±_ Zxm ^ t (2 . 35) 

Their perturbation expansion is again performed by expanding with respect to A, 
(12.3ip . using (I2.33P and applying the Wick rules. Feynman diagrams are the pictorial 
representation of the perturbation expansion for the interacting n-point functions. 
They involve tree and loop diagrams; see Figure Q] for a few examples. A diagram 
is called connected if all outer lines are connected to each other. In the functional 
calculus, the connected diagrams are represented by the generating functional W\(j) = 
\og Z\(j). With this in mind, in what follows we shall always restrict attention to 
connected diagrams. For more details and further reading, we refer for example to [U 
[231 EHl [22]. 
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X-2 




X 4 £3 




X4 X4 X5 Xq 




Xl X\ X\ X 2 X\ X 2 X\ X 2 X3 

(a) (b) (c) (d) (e) 

Figure 1. Feynman diagrams for the interacting ra-point functions. 



For the later comparison with the classical theory, we now give the combinatorics 
of the diagrams. We use the short notation 



A F (j,j)= A F (x-y)j(x)j(y)d 4 xd 4 y. 



(2.36) 



Lemma 2.1. Suppose that we consider a contribution to G\(x\, . . . ,x n ) of order k in 
perturbation theory. Then G\ vanishes unless n is even. Then the contribution can be 
represented by a Feynman diagram involving 



N = - + 2k 
2 

lines. Analytically, the contribution can be written as 

N 



(-i\) k i N f Ai ••• / A fc fl A H---) 



1 

For connected tree diagrams, we have the additional relation 

n 



k 



1 



(2.37) 



(2.38) 



(2.39) 



Proof. According to (l2~35i (l2~3B and flISfy . 

5 5 



G\(xi,...,x n ) = {-i) r 
k 



1 



iX 
"1" 



Sj{xi) 6j(x n ) 
8 n4 



63 (yi 



d yk 



. 8 \4 1 
''Sj(y k )) Nl 



N 



i=o 



For this contribution to be non-zero, we must have as many derivatives 6/Sj as fac- 
tors j. This gives (I2.37p . Carrying out the derivatives, all the factorials and factors 1/2 
are compensated by combinatorial factors. Collecting the factors of —i and using (I2.37|) 
gives QZT5ED . 

To prove (I2.39D we note that in the case k = 0, only the two-point function con- 
tributes. Obviously, each vertex increases n by two. □ 

We finally recall how to take the "classical limit" in which only the tree diagrams 
remain. As we want to take the limit K — > 0, we clearly need to work in more general 
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units where h 7^ 1 (but still c = 1). Denoting the length scale again by I and the mass 
scale by m, all objects have dimensions of £ and m; more precisely, 

[S] = [h] = m£, M = yf, W = ^« ^] = m - ( 2 - 4 °) 

Inserting a factor of into the exponent in (|2.32p and computing the path integral, 
one sees that (|2.33p is to be modified to 



Z [j] = exp I -— J J j(x) A F (x - y) j(y) d 4 xd 4 yj . (2.41) 
Defining the interacting n-point functions by 

G x ( Xl , ...,x n ):= [ <j)(xi) ■ ■ ■ 4>{x n ) e^ s , (2.42) 



a short calculation using (|2.4ip shows that (12.38H is to be replaced by 

N 



H n + Sk-N { _ iX)ki N f ^ f d 4 yk -Q Ap{ 

J J f) 1 



1 

(note that every factor <p gives the functional derivative h5/5j and that expanding the 
exponential in (|2.42p gives powers of X/h). Using (|2.37p . we conclude that 

G x (x 1 ,...,x n )^h k+ ^ . (2.43) 

For connected tree diagrams, we conclude in view of (|2.39p that G\(x±, . . . , x n ) ~ K n ~ l . 
In order to understand this scaling, one must keep in mind that for describing a 
scattering process, one needs to replace the outer lines by the incoming or outgoing 
fields. Truncating the diagrams gives a factor h~ n (because in view of (|2.43p every 
outer line Go(x,y) carries a factor H). Hence the truncated diagrams all scale ~ H . 
After rescaling the diagrams by a common factor h, in the connected tree diagrams all 
factors of h drop out. But (|2.43p also shows that every loop gives rise to an additional 
factor h. Thus taking the limit h — > 0, precisely the tree diagrams remain. 

3. A Nonlinear Classical Scattering Process 

In this section we shall analyze a scattering process for a classical field (ft being a 
solution of the classical field equation (|2.ip . The physical situation which we have in 
mind is shown in Figure [2j We have an incoming wave which may be composed 
of several wave packets. We assume for convenience that for large negative times, the 
wave packets become more spread out in space, implying that their amplitude becomes 
smaller. This has the advantage that in the limit t —> —00, <f> can be treated as a free 
field (making it unnecessary to "switch off" the interaction adiabatically) . Similarly, 
for large positive times, 4> should evolve into a free field <p OU f We assume that the 
interaction takes place inside the strip / = [— T, T] x R 3 for a given large time T (in 
the end, we will take the limit T — > 00). Inside the interaction region, (ft satisfies the 
nonlinear field equation (|2.ip . whereas outside the interaction region, <p is a solution 
of the linear wave equation. 

In the following subsections we first summarize how the Hilbert space of classi- 
cal fields can be regarded as a subspace of the quantum Fock space. Afterwards, 
we discuss the perturbation expansion of the classical (non-linear) Cauchy problem. 
The corresponding pictorial representation of this expansion significantly differs from 
Feynman diagrams obtained in perturbative quantum field theory. To overcome this 
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> 



Figure 2. A nonlinear classical scattering process. 

difference we motivate and introduce the so-called "nonlinear classical measurement 
process" . This will make it possible to introduce the classical scattering operator and 
the classical n-point functions. The pictorial description of the perturbative expansion 
of these classical n— point functions is then compared with the tree-diagrams of the 
corresponding quantum field theory. 

3.1. The Free Classical Field. In preparation, we need to describe free fields, being 
solutions of the linear wave equation 

□0o = . (3.1) 

We assume for simplicity that 0o is smooth and spatially compact (note that, due to 
finite propagation speed, the last property is preserved under the time evolution), 

0o 6 C 00 (R,Co°(R 3 )) . 

For free fields, the coupling constant A vanishes, and thus the classical energy (|2.2[) is 
quadratic. Thus by polarization we can introduce a scalar product, 

(00, ifo) := - / (Wo + V0 O • Wo) d 3 x . (3.2) 

Due to energy conservation, this scalar product is independent of t. It makes the space 
of free solutions to a Hilbert space (Jio, (.,-))• 

It will be useful to represent the function 0o i n various ways. First, in view of 
the unique Cauchy development, we can describe 0o by its initial values at any given 
time t, 

*(t) := (0o, ^0o)|* G C^R 3 ) x C °°(M 3 ) . (3.3) 
Next, we can represent 0o similar to (|2.7p and (|2.8p as a Fourier integral supported on 
the upper and lower mass cone, 

0o(x) = 7^)4 / ^ (Mk) e-*» t+Ck * + Mk)e^- rk£ ) (3.4) 

where again uj = \k\. Then the energy scalar product (|3,2p becomes 

'*>■*>> = wtS^tM^- (3 - 5 > 
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d k / — * — * ~ — \ 

— [Mk) Mk) + <h(k) Mk)) ■ (3.7) 



This scalar product has the same units as the classical energy (see after (|2.2p ). 

[(0o|^o)] =r l - (3.6) 

Another scalar product can be obtained by identifying the classical solutions with 
vectors of the one-particle Fock space: For given <j>o, we seek a Schwartz function / 
such that 

0o(x) = 2Re(O|0 o (x)|/>. 
Using (|2.2U|) and the commutation relations, we find that / is uniquely determined by 

f(k) = Mk) ■ 

We introduce the real scalar product (.|.) on 3t~o by (^o\4>o) := (/)/)• Using (|2.2ip . 
this scalar product has the representation 

1 f d?k ~ 

1 f d 3 k 
~ 2(2^)3 

This scalar product is Lorentz invariant. Moreover, comparing with (|3,5p and ()3.6|) . 
one sees that it is dimensionless, 

[(MM =i°. (3.8) 

It is easy to verify that, up to a multiplicative constant, the scalar product (|3.7p is 
indeed the only Lorentz invariant scalar product which can be introduced on the free 
classical scalar fields. We thus obtain the isometric embedding 

(:Ko,(.|.>W(J,(.|.» 

* \f) • (3 ' 9) 

3.2. Perturbative Solution of the Classical Cauchy Problem. We now return 
to the nonlinear wave equation (12, ip . For simplicity, we only consider solutions which 
are smooth with spatially compactly support, 

4> e C DO (R,Cg°(R 3 )) . 

We denote the set of all such solutions by %. Note that, since our equation is non- 
linear, % is not a vector space. Moreover, there is no bilinear form or scalar product 
on H. The only quantity available is the classical energy (|2.2|) . 

E : U -»■ . 

For the description of the Cauchy problem, it is convenient to again combine <j> and 
its time derivative to a two-component function <I>, (|3,3p . and to write the equation as 
a system of first order in time, 



**=*w *(*)=U0-V/6j- (3 - 10) 

Then the initial data at some time to is a vector 

$(to) G C °°(R 3 ) 2 . 

The time evolution obtained by solving the Cauchy problem gives rise to the mapping 
U(t,t ) : C^°(R 3 ) 2 -»• C^iE 3 ) 2 with <f>(t) = U(t,t ) *$(t ) . 
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The operator U(t,to) is referred to as the time evolution operator. Here the star 
emphasizes that it is a nonlinear operator. But clearly, the time evolution has the 
group property, 

U(t",t')*U(t',t) = U(t",t) for all t, t' , i" G R . (3.11) 

Let us rewrite the solution of the Cauchy problem perturbatively. We decompose 
the operator Sj into its linear and nonlinear parts, 

= j3o$ + A£($) , 

where 

In the case A = without interaction, the equation (]3.10p is linear. The correspond- 
ing Cauchy problem can be solved formally by an exponential, 

$(t) = e (*-*o)J5o _ ( 312 ) 

In the next lemma, we express the linear time evolution operator e(* — in terms of 
the retarded Green's function (j2. 17|) . 

Lemma 3.1. For any t > 0, the operator e ts ^° can be written as 

( e *»°$)(z0 = J R t (x-y) $(y) d 3 y , (3.13) 
where the integral kernel is the distribution 

R tW - {-AS A (t,x) -d t s A (t,x)J- (d - i4j 

Proof. From the representation (|2. 18j) one sees that for positive times, S A (t,x) is a 
solution of the wave equation. Thus a short calculation yields that the function (|3. 13() 
is is a solution of the equation {dt— 5)o)(e tf)0 &) = 0. Next, one easily verifies from f|2. 18j) 
that 

lim S A (t, x)=0 and lim d t S A (t, x) = -5 3 (x) , 

showing that (I3.13P has the correct initial values at t = 0. □ 
In order to treat the interaction perturbatively, we write ()3.10p as 

(d t -Sjo)$ = XB($). (3.15) 
Similar to the interaction picture, we set 

*iBt(t) = e"* 8 " $(t) , BiBt(^mt) = e-** B(e m ° <& int ) . 
Then (13.15P simplifies to 

<9 t $ int = \B int (<f> iat ) • (3-16) 

Making an ansatz as a formal expansion in powers of A, 

$- mt (t) = *W + A*g(t) + A 2 $S(t) + • • • (3.17) 

(where is a time independent wave with the correct initial data), we obtain for n = 
1,2,... the equations 

**Sm = -Gj2 (x) ). < 3 - 18 > 
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Figure 3. The tree diagrams in the perturbation expansion of the 
classical Cauchy problem. 



where we set x = (t,x) and 
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,-tSjo 







p [n \x) 



p( n \x) J ' 

E 



-0( a )(x)0( fe )(x)^(x) . 



a,b,c with a+6+c=n— 1 
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(3.19) 
(3.20) 



Integrating (|3.18p on both sides and using that $j nt (to , s) = for n = 1,2, 
obtain 

Transforming back to $ gives 

Using the explicit form of the free time evolution operator (13.14p . we obtain 
$ {n) (t,x) = - J'dr J d 3 yR t . T (x-y)^ p{n) {T ^^j . 
Setting y = (r, y), the perturbation expansion (|3.17p becomes 

oo 

0(x) = ^AV n) , 



we 



(3.21) 



n=0 



where (j)^ 
and 



is a solution of the free wave equation with the correct initial values, 



0(«)(x) = / S A (x,y) pW(y) d 4 y (n > 1) . (3.22) 
J{y°>to} 

We thus obtain an expansion of in terms of tree diagrams, which can be depicted as 
shown in Figure [3j 

Let us briefly compare the above perturbation expansion for the classical time evo- 
lution with the perturbation expansion for the quantum field as outlined in Section [2T2l 
The equation for the classical evolution equation in the interaction picture (|3.16p can 
be regarded as the analog of the quantum time evolution operator U in (|2,24p . Al- 
though formally similar, they differ in that the classical evolution is non-linear, whereas 
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the quantum evolution is linear. The resulting perturbation expansion for the classical 
field differs from that for the quantum field in that only non-loop diagrams appear 
(see Figures [T] and [3]) . Moreover, the diagrams corresponding to the perturbative ex- 
pansion of the classical field always involve the retarded Green's function, instead of 
the Feynman propagator, and they all have exactly one outgoing leg. 

3.3. A Nonlinear Classical Measurement Process. Let us try to mimic the con- 
struction in linear scattering theory as outlined in Section [2TT1 (see (|2.3p - (|2.6p ). De- 
noting the free classical dynamics f|3. 12|) by Uo(t,t') = e ( - t ~ t ^°, we can introduce in 
analogy to (|2.5p 

n ± iP := U(t , ±T) * U (±T, t ) V : -> U , 

where at time T we identify free and interacting solutions. Note that the operators £l± 
are nonlinear. Moreover, as H is not a linear space, we cannot take their adjoints. 
Rewriting (|2.6|) with the inverse, the obvious idea is to define the nonlinear scattering 
operator by 

(ft-)- 1 * n + = U (t Q , -T) * U(-T, T) * U (T, t ) , 

where in the last step we used the group property (|3. 1 1 j) . As the operator [7o(T, t ) 
maps a free solution at time to to the same solution at time T, it is the identity on "Kq. 
Thus the naive ansatz for the scattering operator is 

S : H -> ^0 , S <t>out = U(-T, T) * out . 

Similar to (|2.28p . a matrix element of the scattering operator would be given by 

(Pout,a in ) = (S- 1 *f3 in ,a in ) (3.23) 

(alternatively, one could work with the scalar product (|3.7|) : this would make no dif- 
ference for the following consideration). However, this naive approach does not work 
for the following reasons. The first problem is that the construction manifestly dis- 
tinguishes a direction of time. This can be seen from the fact that, similar to (|3.22p . 
the perturbation expansion for U(—T, T) will involve only advanced but no retarded 
Green's functions. Moreover, as U(—T, T) is a nonlinear operator, the expression ()3.23j) 
is nonlinear in /3; n , but it is linear in a- m . A related problem is that S is not unitary, 

(S^ 1 -k /3 in , a in ) / (/3 in , S -k a in ) . 

This inequality is obvious because the left side is linear in a m , whereas the right side 
is not. 

In view of these problems, it is not obvious conceptually how to introduce a non- 
linear scattering operator. In order to clarify what we have in mind, we first discuss 
the physical setting in more detail. As explained at the beginning of Section [3l we 
consider a solution <j) of the nonlinear wave equation (|2.ip which as t — > ±oo goes over 
asymptotically to free solutions 4>\ a and cj) out (see Figure [2]). Our wave is classical. 
First of all, this means that the field is not a quantum field in the sense of second 
quantization. Moreover, the usual point of view is that a classical field can be deter- 
mined at every space-time point by measuring its effect on a point-like "test particle" , 
which satisfies the classical equations of motion. By making this test particle suffi- 
ciently small, one can measure the classical field to any precision without changing the 
system (note that this is different in a quantum system, where a measurement changes 
the system by projecting the quantum state to an eigenstate of the observable). In 
particular, one can determine <ft at some initial time to and seek for a solution at later 
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times by solving the Cauchy problem (see Section I2.1|) . The unique solvability of the 
Cauchy problem shows that the classical system is deterministic. Moreover, the fi- 
nite propagation speed shows that causality is respected. In a perturbative treatment, 
causality becomes manifest by the fact that only retarded Green's functions appear 
(see Section [372]) . 

In what follows, we modify the above setting in that we do not want to refer to 
test particles. This is because our model only involves the bosonic field cj), but there 
are no classical point particles which could be used for measurements. Moreover, 
we do not want to make the idealized assumption that the system can be measured 
without disturbing it. More specifically, the only observables which we want to use for 
measurements are differences of classical energies. Working with the classical energy 
seems natural because it is a conserved quantity of the system. We restrict attention to 
energy differences because a constant offset to the classical energy can be interpreted 
as being the contribution of a background field which cannot be detected. With these 
restricted measurements, the observer cannot determine the field at every space-time 
point. In particular, he cannot determine the field <f> pointwise at an initial time to. 
Therefore, he cannot solve the Cauchy problem to make predictions on the future. 
Thus, although our system is still classical and deterministic, the observer cannot 
predict the outcome of an experiment with certainty. 

We now discuss how a scattering process can be described in this restricted frame- 
work where only differences of classical energies are allowed for measurements. Suppose 
that an observer at some large time T wants to detect the result of the scattering pro- 
cess. The most obvious method to make observations is to modify the wave cj) by some 
"test wave" 5cj> and to consider how the classical energy 112. 2|) changes. If the test wave 
is taken into account linearly, the energy is perturbed by 



Using that at time T, the outgoing wave packets should be so spread out that the term 
involving A can be dropped, we obtain the simple expression 



where (., .) is the scalar product obtained by polarizing the free energy (|3.2j) at time T. 
Despite the fact that we are considering purely classical fields, the resulting situation 
has similarity with a measurement process in quantum mechanics. Namely, the com- 
putation of the expectation value in (I3.24p can be interpreted that a measurement 
of the wave cj) is performed with a prepared "end-state" 5(f). By modifying 5<j), one 
can determine cj> completely. In particular, one can measure the distribution of cj) in 
momentum space. 

The important point for what follows is that the relation (]3.24p only holds in the 
linear approximation. If the amplitude of 5(j> is increased, it has a nonlinear effect 
on cj>, which influences the result of the measurement process. In order to analyze 
such nonlinear effects, we consider the situation shown in Figure [H It is convenient 
to describe the incoming field by an inhomogeneity p- m which lies in the past of the 
interaction region and generates a wave moving to the future. Similarly, the observer 
is described by an inhomogeneity p ut> which is located in the future of the interaction 
region and generates a wave moving to the past. The resulting wave cj> will be a solution 





(3.24) 
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Figure 4. A nonlinear classical measurement process. 

of the nonlinear inhomogeneous wave equation 

_ ncf> - -<f> 3 = \ ( Pin + Pout ) . (3.25) 
o 2 

Due to the inhomogeneities, the classical energy (12. 2p is not conserved. But it is 
conserved in the regions of space-time where both p- m and p out vanish. We can thus 
distinguish three different energies: the energy E in the interaction region, the en- 
ergy E_ in the past of /?i n , and the energy E + in the future of p out . 

As explained above, we take the point of view that the observer can only measure 
energy differences. Moreover, as we are interested in the effect of the interaction, it 
is convenient to also subtract the corresponding energies of the free fields, which we 
denote by the additional superscript "free" . This leads us to introduce the quantity 

AE := (E - E iree ) - (E + - E% ee ) . (3.26) 

It gives the energy shift caused by the interaction as being detected by the observer 
at late times. We interpret AE as the result of our classical measurement process. 

3.4. The Global Perturbation Expansion. In what follows, similar to (|2.36p we 

use the short notations 

S A (x,p)^ J S A (x,y)p(y)d A y, S A ( Pl ,p 2 ) = J J S A (x,y) p x (x) p 2 (y) d A x d A y , 

also analogously for other two-point distributions. In order to model the classical 
measurement process (see (|3.25p and Figure H]), we choose the free solution </>(°) as 



= 0in + 4>out With 

1 1 (3-27) 

0in(aO := - S A (x,p in ) , 4>out{x) := - S v (x,p out ) . 

It is a solution of the inhomogeneous linear equation — D^o = (pin + Pout)/2- In order 
to construct a corresponding solution of the nonlinear equation (|3,25p . we make a 
perturbation ansatz formal power series in A, 

oo 

(f> = A n {n) • (3.28) 

n=0 
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By direct calculation one verifies that the functions <jP~\ 4>( 2 \ ■ ■ ■ can be written itera- 
tively as 



= J S^(x,y) p^\y) d% with 

1 (3.29) 

P (n) (y)--=xi £ ^ (a ¥V c) )(y), 

a,b,c with a+b+c=n— 1 

where S^ 71 ' is a general Green's function, meaning that it is a distributional solution 
of the equation 

-DS^(x) = 5 4 (x) . 

Here xi 1S the characteristic function defined by xi{ x ) = 1 if x € J and xi{ x ) = 
otherwise. It has the effect that in (|3.29p we only integrate over the interaction region. 
This is a technical simplification which will be used in the proof of Lemma 13.31 It 
has no physical significance because in the scattering process under consideration, 
we assumed that the dynamics outside the interaction region is linear, so that cubic 
expressions in can be omitted. 

Since we want to be a real-valued function, the Green's function should be real. 
This gives rise to the general ansatz 

SW(x, y) = S (x, y) + c n P (x, y) + id n K (x, y) , (3.30) 

where Sq is the causal Green's function 

So(x, y) = I (S^ + S*) (x, y) = f ^ e -*(-») , (3.31) 

whereas Kq,Pq are the fundamental solutions of the free wave equation (|2.15|) and 

^ 6(k>) , (3.32) 

and the constants c n and d n are free real parameters (note that Po is real-valued, 
whereas Kq is purely imaginary). Thus the perturbation expansion involves two real 
parameters, which can be chosen freely to every order in perturbation theory. 

At this stage, there seems no reason why the parameters c n and d n should depend 
on the order n of the perturbation expansion (note that this is different for the pertur- 
bation expansion for fermions, where the concept of the Dirac sea and normalization 
conditions give rise to a non-trivial combinatorics of the operator products OH]). 
Therefore, in the remainder of this paper we always make the simplest choice where 
all parameters c n and d n vanish, 

S^(x,y) = S (x,y). (3.33) 

The resulting perturbation expansion (|3.29p and (|3.33p is neither advanced nor re- 
tarded and should thus be regarded as a global expansion in space-time. Therefore, 
we refer to it as the global perturbation expansion. It gives rise to a nonlinear map- 
ping p which to every associates a corresponding solution <f> of the nonlinear 
equation, 

p : ?{ (0) -> U . (3.34) 

We close this section with a few remarks. We first point out that our choice of 
Green's functions (|3.33p was mainly a matter of simplicity and convenience. But more 
general choices of the coefficients c n and d n are also possible and might be worth 



20 



F. FINSTER AND J. TOLKSDORF 



considering in the future. We also remark that the vanishing of the coefficients d n 
in (|3.30p can be motivated by time reflection symmetry as follows: Note that in our 
classical scattering process we did not distinguish a direction of time, meaning that 
the setting is symmetric under the transformations 

t — > — t and 0j n ou t (3.35) 

(clearly, this transformation exchanges the observer in the future with an observer in 
the past, and also leads to the replacement E + -R- EJ). Therefore, it seems natural to 
impose that also the perturbation expansion should be invariant under ()3.35p . Since 
the distributions So and Pq are symmetric in its two arguments, whereas Kq is anti- 
symmetric, this leads us to impose that d n = 0. 

We finally note that the Feynman propagator (j2.30p can be written as 

Ap = So — iirPo . 

This is of the form f|3.30j) . but with the parameter c n = —m being imaginary. This 
shows once more that the Feynman propagator is complex-valued (indeed, it is imagi- 
nary in spacelike directions). Therefore, it cannot be used for the perturbation expan- 
sion of real- valued classical fields. 

3.5. A Perturbation Expansion of the Classical Energy. We now consider the 
classical energy (|2.2p of our solution <j>, 

E (t) = I (U' + l IWI 2 + ^ </>*) d'x • (3-36) 



Jt=const \2 2 4! j 

As in Figure HI we denote the energy in the interaction region by E, whereas E± are 
the energies for large positive and negative times, respectively. 

Proposition 3.2. The classical energies E and E± have th> 

m ■ , \ A f 

E+ = —K {p in + p, p in + p) + hm — / 

4 t-+oo 4! J t=const 

_ m r > / „ _n A 

E = —K [p in - p out + p, p in - p out + p) + — 

= -jAo^n - p ou t ~ p, Pin ~ Pout ~ P) + 

m ■ , s A 

E- = —Kolpout + P, Pout + P) + hm — 

4 4! Ji=const 

where Kq(x,u) = d x oKo(x,y) is the time derivative of the distribution (j2. 15[) . Here 
the function 4> is defined inductively by ()3.27p ()3.29[) . and p is given by 

oo 

p = J> (n) - ( 3 - 41 ) 

n=l 

Proof Due to energy conservation, the energy E can be computed at any time t G 
[— T, T]. We first compute it at time t = T to derive (|3.38p . We again polarize the free 
energy to obtain the scalar product (.,.) introduced in (|3.2p . We now substitute the 
perturbation expansion (|3.28|) into (|3.36p . Multiplying out the free part of the energy, 
we need to compute expressions of the form 

(S v (.,x), S (.,y)) , (5 A (.,x), 5 (.,y)) , (5 (.,x), S (.,y)) 



2 expansions 




t> 4 d 3 x 


(3.37) 


[ 4> 4 d z x 


(3.38) 


Jt=T 




[ 4 d 3 x 


(3.39) 


Jt=-T 




4> 4 d 3 x , 


(3.40) 
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Furthermore, we know that the argument x of the factors S A (.,x) and Sq(.,x) always 
lies to the past of time T (because x £ I U supp/9i n ), whereas the argument of the 
factors S v (.,x) lies to the future of the time T (because in this case, x G supppout))- 
According to (|2.16p . f|2. 18[) . (|3.3ip and (|2.15p . we may thus exchange the Green's 
functions with the replacement rules 

S A (.,x) ->■ -2mK (.,x) , S v (.,x) 2m K Q (.,x) , S (.,x) ->■ -m K (.,x) (3.42) 

by corresponding fundamental solutions Kq. After these replacements, we can apply 
Lemma 13.31 below. In particular, we obtain for any n,n' > 1, 

= (_^)2 | j (K (.,x), K (.,y)) pW(x) p«\y) d 4 xd A y 

= -* 2 J [ K,{x,y) ,<»'>(y) d 4 xd 4 y = | ^(pM ,<»'>) . 

Adding all the terms of the perturbation expansion gives ()3.38p . 

To derive the expansion (|3.39p . we compute the energy similarly at time t = —T. 
Then the interaction region lies in the future of t, and therefore the replacement rule 
for So in (I3.42|) is to be modified to 

S {.,x) -> iirK (.,x) . 

The expansion for E + is derived similar to f|3.38j) . However, as the argument x of 
the factors S v (.,x) now lies in the past of t, all the term involving a scalar product 
with S v (.,x) drop out. Then energy conservation allows us to take the limit t — > 
oo. Similarly, the expansion for E- is obtained from (|3.39p by omitting the terms 
involving S A (.,x) and taking the limit t — > — oo. □ 

It remains compute the energy scalar product of two factors Kq. 
Lemma 3.3. The following relation holds in the sense of distributions, 

(K {.,x), K (.,y))=~k (y,x). 
Proof We write the scalar product in momentum space as 

(K (.,x), K (.,y))=lim \ j ' d 3 z j ^* e(k°) S(k 2 ) e' ik °^ 

x J e(g 0) gtf) e - iq °(r-y°) + i<?(z-g) (^_ k q _ ^ 

Carrying out the spatial integral, we only get a contribution if q = —k. Setting Q = q°, 
we thus obtain 

1 f d 4 k f°° 

(K (.,x),K (.,y)) = l^-j-4sj dne(k°)8(k 2 )e(n)S(n 2 -\k\ 2 ) 

Due to the two delta distributions, the integrand vanishes unless Q = ±k°. In the 
case f2 = the last factor vanishes, because 

-k°n + \k\ 2 = (k°) 2 -\k\ 2 = k 2 = o. 
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Thus we may set = —k° to obtain 

1 f d^k i 



(K (.,x), K (.,y)) =-±j ^ ^ i(^) e-^^H^GMO 2 |£| 2 

Comparing with (|2.15p gives the result. □ 

We now analyze the result of Proposition [321 First, if we assume that the solution (j> 
is dissipative and for large times goes over to a solution of the free wave equation, the 
last integrals in (13.37P and (]3.40j) vanish. This dissipation could be made precise 
using the techniques described in [THJ EJ [TTl EQ] . Here we take the decay property for 
granted and simply drop the last integrals in (13. 37ft and (|3.40p . Similarly, if we choose 



the interaction time T sufficiently large, we can also drop the last integrals in (13.38| ) 
and (|3.39p . Then taking the difference of ()3.38p and (|3.39p gives the identity 

K (p in - Pout, P~) = . (3.43) 

Using this identity, we obtain 

ITT 

E+ = —K [p in + p, p in + p) 

E = —K (p in - p out , p in - p out ) + — K (p, f>) 
E- = —K (pout + P, Pout + P) ■ 

Here the function p describes the component of the wave generated by the interaction, 
whereas p- m and p Q ut describe the free incoming and outgoing components. In order to 
analyze how the energies are effected by the interaction, we subtract the free energies 
as given by 

^ ee = ^^ (pin,Pin) 

E bee = — ii"o(pin — Pout, Pin ~ Pout) 



El ee = — -Ko(Pout,Pout)- 

We thus obtain the compact formulas 



E-E^ = ^k (~p,p) 



E + - E%™ = E„- = y k (pout , p) + ^K (p, p) . 



Let us discuss these formulas. First, the formulas for the free energies can be 
understood immediately from the fact that for large positive (negative) times, only 
the wave generated by p; n (resp. p ou t) is present, whereas in the interaction region, 
both waves are superimposed. The interaction affects the energies by two different 
contributions: The term 

Z -jK (p,p) 

is delocalized in space-time in the sense that it adds to the energies E± and E in exactly 
the same way. As a consequence, this contribution is not accessible to the observer 
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Pin Pin Pin Pin Pin Pin Pin Pin Pin Pin Pin Pin 

Figure 5. The perturbation expansion of the energy shift AE in the 
classical measurement process. 

who can only measure the energy difference E — E + , whereas a constant offset to the 
energy is not detectable. The situation is different for the contribution 

yi^o(pout,p), 

which only effects the energy in the interaction region. The observer at time T can 
detect this contribution by measuring the classical field energy before and after modi- 
fying the field by p ou t- We now identify this contribution with (|3.26p and bring it into 
a more convenient form. 

Proposition 3.4. In the classical measurement process described in Section \3. 31 the 
energy shift (|3,26p as measured by the observer at time T is given by 

AE = \ S v (p,p out ) = -\ S A (p, Pin ) . (3.44) 

Proof. Comparing the above formulas for the classical energies with (|3.26j) . we find 

AE = yiTo(pout,p) = y^o(pin,p)- 

The function p is localized in the interaction region. Therefore, we can use (|3.42j) to 
replace K by the time derivatives of the causal Green's functions. □ 

We remark that (13.440 shows that our constructions are symmetric under time reversals 
(note that a time reversal corresponds to the replacements p ou t ^ Pim S v <-> 5 A 
and d t —d t ). 

Expressing p by (13.41 b and the perturbation expansion (|3.27p - (|3.29p . we obtain an 
expansion of AE in terms of tree diagrams. The resulting diagrams are shown in 
Figure [5j This perturbation expansion differs from that for the solution of the Cauchy 
problem (see Figure [3j) in several ways: First, we can now have several outgoing 
legs. Second, the inner lines now involve the Green's function So (instead of the 
retarded Green's function 5 A ); this is a direct consequence of our choice of the Green's 
function (I3.33P in the global perturbation expansion. Finally, the time derivative 
in (|3.44p has the effect that in each diagram, exactly one of the outgoing legs is 
differentiated with respect to t. This leg is denoted by a dotted line. This time 
derivative obviously destroys the Lorentz invariance. This can be understood from the 
fact that as an energy, AE is not a Lorentz scalar but the zero component of a Lorentz 
vector. 
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3.6. The Classical Scattering Operator and Classical n-Point Functions. Let 

us reconsider the previous constructions. The result of the classical measurement pro- 
cess was associated to a suitable difference of classical energies A.E as given by (|3.26[) . 
To compute these energies, we started from the free solution (|3.27p . performed the 
global perturbation expansion, and could then evaluate the energy integrals (see Propo- 
sition [3T2]) . Describing the perturbation expansion by the nonlinear operation (|3.34[) . 
we thus obtain the nonlinear mapping 

AE : IK x 5C -> R , AE(^ out , </> in ) := AE{p * (0 out + in )) . (3.45) 

Plugging in the perturbation expansion (|3.29p and multiplying out, we get an infinite 
sum of multilinear mappings, 

00 

AE(4> outl (p- m ) = AE p!q ( 4> out , 4> ou t , </> in , . . . , 0m ) ■ (3.46) 
— n S v ' v y ' 

P'1~ u p arguments q arguments 

At this point, it is helpful to apply the standard procedure of rewriting the multi- 
linear operators as linear operators on the tensor product. To this end, we introduce 
the free p-particle classic bosonic Fock space 

^dass = ffio g • • • ^oj symm , 
n factors 

where "symm" means total symmetrization (and 3® are simply the real numbers). 
On 3^i ass we let (.|.) be the scalar product induced from the Lorentz invariant scalar 
product (|3.7j) on "Kq. Taking the direct sum gives the 

00 

classical Fock space 3" c iass = -^dass • 

n=0 

We remark that lifting the isometric embedding (|3.9p to the tensor product yields a 
canonical isometric embedding of the Fock spaces, 

(J'class, (•!•)) ^(?, (•!•)) • 

The two Fock spaces differ by the fact that Sciass is a real, whereas 3 is a complex 
vector space. In fact, 3~ is the complexification of the image of 3~ c iass under the above 
embedding. 

As the multilinear mappings in (|3.46p give rise to unique linear mappings on the 
corresponding symmetric tensor products, we obtain a unique bilinear operator 

AE : J dass x J class -> R . (3.47) 

The classical scattering operator S is introduced by 

S : J dass -> J class with AE($ ollt , * in ) = -(d t $out I & $ in ) • (3.48) 

In the next lemma we bring the classical scattering operator into a more explicit form. 
Our proof will also explain why the factor dt in (|3.48p is needed in order to make the 
definition Lorentz invariant. 

Proposition 3.5. The scattering operator has the matrix elements 

(4>L ® • • • ® font I ^ <tL ® ■ ■ ■ ® O = W^out, • • • , Ct^in, ■ ■ ■ , O • (3-49) 

Here the functions S p + q have a perturbation expansion in terms of Feynman tree dia- 
grams. The order k and the total number of lines N are given by (|2.37p and (|2,39p . 
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where n = p + q. The contributions to S n without the outer lines are given analytically 
by 

/,. N—n 
Ai ••• / n ( 3 - s °) 

(where So is the Green's function (|3.3ip ), 

The functions S n (xi, . . . ,x n ) are referred to as the classical n-point functions. 
Proof of Proposition \3.5[ Combining (|3.44|) with (|3.27p , we obtain 

AE = \J p(x) <p on t(x) d 4 x . (3.51) 

Expressing p by the perturbation series (|3.4ip and (|3.29[) . one gets an expansion in 
terms of tree diagrams. This explains (|2.37p and ()2.39p . In order to get hold of the 
combinatorics, we consider the summand AE p q in (|3.46p and number the arguments 
by ^outi ■ • • > ^out ano - ^in' ■ ■ ■ ' ^fn- Moreover, we restrict attention to the contributions 
to AE p ,q corresponding to a fixed tree diagram, whose outer lines are labeled by the 
outer arguments. According to (13.51 f) . exactly one of the outer arguments 4>\ n , . . . , 0f n 
involves a time derivative. If we fix this dotted outer leg, the only combinatorial factor 
comes from the freedom to permute the three incoming legs at each vertex. This leads 
to a factor 3! at each vertex, which cancels the factor 1/6 in (|3.29p . We thus obtain 

= 2 ( 5 'n( ( Aout> • • • > dt> ■ ■ ■ ) 40 + H Sn(4>luti • • • ) <^ut> • • • > </>iS) 

= \ (WLt ® • • • ® €ut) I ^ <i>L ® • • • ® O > (3-52) 

where S n is formed according to (]3.50p , and £H is given by (|3.49p . Comparing with (j3.48j) 
shows that indeed agrees with the classical scattering operator. This concludes the 
proof. □ 

Let us compare the classical n-point functions with the n-point functions of quantum 
field theory as computed in Lemma 12.11 One obvious difference is that the classical 
n-point functions involve only tree diagrams, whereas the loop diagrams of quantum 
field theory are missing. Before comparing the combinatorics (13.50P with (12.38f) . one 
should note that in (13.491) we inserted free fields 4>f n and (jr 0XLt as arguments. This is 
why in (13.501) there are no propagators for the outer lines (this corresponds to the usual 
procedure of considering the "truncated" diagrams). This difference also accounts for 
the additional factors of i in (12.38[) . Namely, according to (12.371) and (12.39p . for tree 
diagrams the identity n = N — k — 1 holds, and thus the factor i N ~ k in (|2.38p can 
be generated in the classical scattering operator if to every outer line we associate a 
factor of i. Then the combinatorial factors in (13.50P coincide with those in (I2.38p . up 
to an irrelevant prefactor which could be absorbed into the definition of the classical 
scattering operator (I3.48P (the reason for the convention (]3.48p is that we want the 
classical scattering operator to be real). We conclude that on the level of tree diagrams, 
the classical and the quantum scattering agree, including the combinatorics of the 
expansion. The only difference is that in the classical n-point functions, the inner lines 
are described by the Green's function So, whereas in the quantum scattering process 
the Feynman propagator Ap appears. This difference can be understood from the fact 
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that the classical field is real, and therefore its dynamics cannot be described with the 
complex- valued Feynman propagator (see the end of Section 13. 4[) . 

We note that the classical scattering operator also differs from the quantum scat- 
tering operator in that it only involves the connected diagrams. As a consequence, the 
classical scattering operator S : 3~ c i ass — > 3~ c iass is not unitary (this is obvious already 
without interaction, because in the case A = 0, the operator S is the identity on the 
one-particle subspace S^lass but vanishes on the many-particle subspaces). A possible 
strategy to obtain a unitary classical scattering operator would be to take the sum of 
tensor products of S, like for example the expression 

x i 

V— s® ■■■ ® s . 

^ n! v v ' 

n— n f ac tors 

This would give a unitary operator in the case A = 0. However, it is not clear whether 
this operator would also be unitary in the case with interaction. If not, one could take 
its polar decomposition and absorb the non-unitary part into the scalar product of the 
classical Fock spaces. Entering the details of this analysis or discussing its potential 
physical significance goes beyond the scope of the present paper. 

We finally remark how the identity (|3.43p can be understood in terms of Feynman 
diagrams. In the proof of Proposition 13.51 we computed AE with the first equation 
in (|3.44p . The identity (|3,43p makes it possible to compute AE alternatively by the 
second equation in (|3.44p . In analogy to (|3.5ip and (|3,52p . we thus obtain 

AE = ~ j p[x) in (x) d A x (3.53) 

AE M = -\ (<^ ut ® • • • ® Ct I ^ d t (cPl ® • • • ® C)> (3-54) 

Subtracting the last equation from (|3.52p . we find that (I3.43P is equivalent to the 
identity for the S'-matrix elements 

{d t $ out | 5 f $ in > + <$ out | & dt&ta) = . 

This means that the value of the Feynman diagram remains unchanged if both the 
incoming and the outgoing waves are translated in time. It corresponds to the usual 
concept of conservation of quantum mechanical energy in a Feynman diagram, stating 
that the incoming energy, defined as the sum of the frequencies of the incoming waves, 
coincides with the outgoing energy. 



4. Nonlinear Classical Scattering in a Stochastic Background Field 

4.1. A Free Stochastic Field. We begin by analyzing a free stochastic field £. By 
"free" we mean that £ € "Kq is a solution of the linear wave equation. Similar to (13. 3p . 
we can characterize £ by its initial values for example at some fixed time to, 

S:=(f,f)|to6CS°(R 3 )xCo°(R 3 )- 

"Stochastic" means that we prescribe a probability distribution for £ which we assume 
to be Gaussian. We write the corresponding probability measure £)H on the initial 
data as 

®~ = ^e-^ F ^d£Ad£, (4.1) 



BOSONIC LOOP DIAGRAMS IN CLASSICAL FIELD THEORY 
where Z is a normalization constant, and F is a positive bilinear form, i.e. 



27 



F(E,E) = j (E(x),FE(x)) r2 S 



x 



with F being a linear operator on the initial data. Next, it is natural to assume 
that our probability measure is translation invariant, meaning that the operator F is 
diagonal in momentum space. Thus setting 

ri 3 Ic 

1(k) e ikx , (4.2) 



E(x) 

the operator F can be written as 



(27T)* 



(FE)(k) = F(k)E(k), 

where F(k) is a 2 x 2-matrix. The fact that 5 is real-valued and that F acts on a 
space of real- valued functions implies that 



E(-Jfe) = E(Jfe) and F(-k) = F{k) , (4.3) 

where the bar denotes complex conjugation of each vector or matrix component. Using 
Plancherel's theorem, we can write the bilinear form F as 

/d 3 k -, - 

— z {E(k),F{k)E{k)) c2 . (4.4) 

A Gaussian probability measure is completely characterized by the covariance C 
defined as the expectation value of a quadratic polynomial, i.e. in our setting 



C(x,y) 



ay) & 



(4.5) 



The expectation values of higher order polynomials can then be computed using Wick's 
theorem (similar as outlined in Section 12.21 in the context of path integrals or see for 
example [7]). 



Lemma 4.1. For any space-time points x,y E I, the covariance 

d 4 k 



5j) is given by 



C(x,y) 



(2^) 4 
e 



5(e 



J!_ e ~ik(x-y) 
M 



iui(x°-y ) I I -I'M 



,F(k) 



-i 



-ioj 
1 



Jui(x°+y () +2t) 



-IUJ 



I r F{k)-\ | 



tu: 



Proof. In momentum space, the Gaussian measure (14. ip becomes 



DE = — exp 
Z 



d 3 k 



{E(k),F(k)E(k)) C 2)dE 1 (k)AdE 2 (k) 



(4.6) 
(4.7) 

(4.8) 



(2tt) 3 

By the standard calculation rules for Gaussian measures (see for example 0), we 
obtain 

J E(k) ® E(k')^ DE = {2tt) 3 8 3 {k - k') F(fc)- 1 . (4.9) 
Next, we know from (|3.13p that 

e(x)= / (-S (x,(t,y)) -S (x,(t,y)))E(y)d 3 y 
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and thus 



m = [ 



cfizi d 3 z 2 



S°(*,(t,*i))\ ^ (g ) (-S°(y,(t,z 2 )) 



-S"(y,(t,z 2 ))J /& 



Similar to (13,421) . we may apply the replacement rule 

s o( x >y) -> -2TriK (x,y) . 

As in Lemma [231 for the momenta we use the notation k = (tJ,k) and q = (O, g). We 
then obtain 



d A k 

W) 4 



e ikx e(uj) 5{k z ) 



d A q 



x e 



i(u-0)t / ( l0J 



(2vr) 4 
(k) <8> H(g) f 



e- iqy e(O) % 2 



in 
-l 



Now we can carry out the Gaussian integral with fj4.9f> . 



/ 



£(s) = 4vr 2 



d 4 k 
(2vr) 4 



e ifca; e(w) 5(A; 2 ) 



d& p ~iQy°+iky 

2vr 



e e(Sl)6(n-\k\' 



7T 



(27T) 4 



e ikx e(w) 



x _ e ~iO.y +iky e i(u-Q)t 



™),F(k)- 



in 



Q = ±OJ 



Rearranging this expression gives the result. 



□ 



Let us discuss the result of this lemma. We first point out that the covariance is 
obviously real-valued and symmetric under permutations and x and y. Moreover, it 
satisfies the free wave equation, 

C(x,y) = C(y,x) and d x C(x,y) = 0. 

The last equation is obvious from (|4.5p because the factors and are solutions of 
the wave equation. Alternatively, it can also be verified in the expression of Lemma [4. II 
by applying the wave operator and using the factor 5(k 2 ) in the integrand on the right. 
In particular, it is impossible to arrange that C(x,y) is a Green's function. 

In order to obtain a covariant theory, we would like the covariance C to be Lorentz 
invariant. This raises the questions for which choices of F the covariance is Lorentz 
invariant, and which Lorentz invariant distributions C can be realized. We first explain 
what the naive answer to these questions is, and then prove that this naive guess is 
indeed correct. The general Lorentz invariant bi-solution of the scalar wave equation 
is a linear combination of the distributions Po and Kq introduced in (]3.32p and (|2.15p . 
As Ko(x,y) is anti-symmetric under permutations of x and y and Pq(x, y) is real- 
valued, we conclude that C should be a real multiple of Po(x,y). 
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Proposition 4.2. Assume that the covariance (|4.5p of the Gaussian measure (|4.ip is 
Lorentz invariant. Then there is a parameter j3 > such that 

C(x,y) = ^P (x,y) (4.10) 

F(E,E) = ^m, (4.H) 
where denotes the scalar product (|3.7p . 

Proof. Lorentz invariance clearly implies that C(x,y) must be rotationally symmetric, 
meaning that -F(fc) must depend only on Combining this fact with the second 
equation in (|4.3p . we conclude that F(k) must have real matrix entries. Thus its 
inverse can be written as 

*■=(; ')■ ^> 

with real functions a, b, c depending on \ k\. Next, Lorentz invariance implies that C(x, y) 
should be independent of F, meaning that the contribution (|4.7p must vanish, i.e. 



TV®- 1 (?)>„=» f - 



Combining this condition with (|4.12p we conclude that 6 = and that c = u 2 a. Then 
the expectation value in (|4.6p becomes 

Substituting these results into the formula of Lemma 14.11 we obtain 

/d^k 
__ 6ik 2 } e ik(*-y) |w| 

This distribution is Lorentz covariant only if a(|/c|) = (2-/r/3 |/c|) _1 for a suitable real 
parameter j3. This gives (|4.10p . Moreover, 

F(k)~ 1 = — ^r(l V) and thus F(k) = 2vr/3 f 1 ^ JL,). (4.13) 
V ' 2Trf3\k\ V° " ) V0 l/\k\J ^ ' 

Next, we write £ similar to (|3.4p as a Fourier integral over the upper and lower mass 
cone, i.e. 

1 F d^h / -> - — — 

i{x) = (2^F J ^7 \^ h) e ^ +tkX + m eiiJt ~ ikX 

Taking the time derivative, evaluating for simplicity at t = and comparing with (j4.2|) 
we find that 

Si(fc) = (tffc) and E 2 (k) = -± 

Hence the functional F as given by (|4.4p and (I4.13P becomes 

F(= S ) = A/« 1,^,.. (4 . 14) 



2vr 7 (2vr) 3 w 

Comparing with (|3.7p gives the result. □ 
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Recall from (|3.8|) that the scalar product is dimensionless. Moreover, the func- 
tional F, being the argument of the exponential in (|4.ip . must clearly be dimensionless. 
Hence the parameter /3 in (|4.1ip is also dimensionless, 

[P}=£°. (4.15) 

This parameter describes the strength of stochastic background field. In the limit (3 — > 
oo, the covariance (|4.10p vanishes, so that no stochastic field is present. If j3 is de- 
creased, the amplitude of the stochastic field gets larger. 

4.2. The Perturbation Expansion and Classical Scattering. We now want to 
include the stochastic background field in the interaction. To this end, we first add 
the free stochastic field £ to the free in- and outgoing fields in ()3.27p . 

fa = fox + 0out + £ • (4.16) 

The corresponding interacting field is then introduced exactly as in Section \3. 41 by the 
perturbation series (|3.28|) and (|3,29p . Expanding the classical energy exactly as in 
Section |3.5| we obtain in analogy to (|3.45p the nonlinear mapping 

AE : M x}{ x}{„^R, AE(^ out , </> in , £) := AE(p * (<£ out + </> in + £)) . 

After expanding similar to (|3.46|) in a series of multilinear mappings, 

CO 

AE = AEp,q,r(<f>out, ■ ■ ■ ,^out,0in, • • • Axxx, £,•••,£) , (4.17) 



P'1' r ^ p arguments q arguments r arguments 

we carry out the Gaussian integral with the probability distribution (|4.ip . 

A£(^ out ,0 in )= AE P,Mut,<f>L) with 



p,q=0 
oo „ 

AE p , q {cg^K) :=£ / AE P!q>r (^ out ^l,e)dE. 

Now the classical scattering operator is introduced similar as in Section [3.61 by (|3.47p 
and (J33SJ). 

We now outline a general procedure to express the classical scattering operator in 
terms of Feynman diagrams. In order to get hold on the combinatorics, it is preferable 
to first derive the expansion for fixed £, and to integrate over 3 afterwards. Similar 
to (|3.47p . we rewrite the multilinear operators in (|4.17p as linear operators on classical 
Fock spaces, 

AE : 3" class x 3" class x 9~ class ^ ^ j 

A£( ^ out ® • • • ® jout , <fen ® • ' • ® </>in , £ ® • • • ® Q = A ^P,g,r(0out^fn» D ■ 
p factors </ factors r factors 

For a fixed stochastic field H st0 ch G 3", we now introduce the scattering operator £[H stoc h] 
in analogy to (|3.48p by 

^[Sstoch] : 3"class — > 3"class with A£($ out ,$i n ,H stocll ) = -(5 t $ out | 5[S stoch ]^ $ m ) . 
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So Po 



(a) 



(b) 



(c) 



(d) 



(e) 



Figure 6. Examples of loop diagrams for the classical measurement 
process in a stochastic background field. 

We next consider the matrix elements of S[H stoc h] for product states. Following the 
procedure in Proposition 13.51 we obtain 



where S p + q + r are again the classical n-point functions. We thus obtain an expansion 
in terms of tree diagrams, but now with additional external lines labeled by £i, . . . , £ r . 
The analytic expressions for the Feynman diagrams including the combinatorial factors 
are again given by (|3.50p . 

In order to carry out the integral over 3, we first replace the factors £i, . . . , £ r by our 
stochastic field £ to obtain an r-multilinear function in ^. Then the H-integration gives 
a Gaussian integral, which can be carried out with the Wick rules (see (j2.34|) ). We 
thus obtain the sum over all pairings between the factors £, where every pair ^(x) £(y) 
is replaced by the covariance C(x,y) (see (14.51) ). In view of (14.100 . every pairing 
generates an additional line described by the distribution Po(x,y)/f3 (where Po is the 
fundamental solution (I3.32p ), We thus obtain loop diagrams as shown in Figure [U 

4.3. Comparison of Classical and Quantum Loop Diagrams. Let us compare 
this expansion to the perturbation expansion of quantum field theory as shown in 
Figure HJ First of all, the expansions are very similar in that they both involve tree 
as well as loop diagrams. The only difference between the diagrams is that the lines 
of the quantum diagrams are formed of the Feynman propagator, whereas the lines in 
the classical diagrams are formed of either the Green's function So or the fundamental 
solution Po; with a specific combinatorics. 

Let us explain the significance of the parameter (3. We already noted in (I4.15|) 
that it is a dimensionless parameter. In the limit /3 — ¥ oo, the covariance and thus 
the stochastic background field vanish. This corresponds to the classical limit where 
no loop diagrams are present. As a diagram with I loops involves a factor /3~ l , by 
decreasing (3 we can make the loop diagrams larger, giving a more "quantum-like 
behavior." Thus the parameter plays a similar role as Planck's constant h in 
quantum theory. In order to make this connection precise, we again need to consider 
the more general units (|2.40p where h ^ 1. Comparing (|3.2p with (|3.7p . it follows 
that the scalar product has dimensions ml. Thus in order for the functional F 



(0out ® • • • ® €ut I ® • • • ® CrY 4>\n « 

= Sp+q+r (0out> ■ • • ) <^uti • • • ' 0jL> & 
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in (|4.1ip to be dimensionless, we need to choose 

const 



h 

with a numerical constant. In this way, the free parameter /3 _1 which describes the 
strength of the stochastic background field can be related to the parameter h in quan- 
tum theory. 

From the analytic point of view, replacing the Feynman propagator A F by So 
or Pq has no influence on the qualitative behavior of the loops. In particular, the 
classical loop diagrams have the same divergent behavior as the corresponding quantum 
diagrams. The following lemma even gives a quantitative connection between the 
classical and the quantum loops. 

Lemma 4.3. Assume that the momenta k\, . . . , k r with r > 1 are all different. Then 
for almost all momenta q, 

/oo poo r 

dq° A F (h + g) ••■ A F (k r + q) = -in / dq° ^P (fc o + q) ]J S {k b + q) . 
-°° J —OO „_1 i /„ 

a—1 b^a 

Proof. We set u = q°, u) a = and p a = \k a + q\ > 0. Then, according to (|2.30|) . 

A F (k a + q) = lim ■ T — (4.18) 

e\0 (LO + UJ a ) 2 -pl + lS 

= i im J_ ( i i ^ . (4.19) 

e\0 2p a \ {UJ + U a ) -Pa + l£ + ^a) + p a ~ IE J 

We take the product of r such factors for fixed e > 0. Since the momenta fej, . . . , k r 
are pairwise different, the resulting function has simple poles for almost all q. In the 
case r > 2, we can close the contour either in the upper or in the lower half plane. 
Computing the mean value of these two contour integrals and taking the limit e\0, 
we obtain 



/oo r 1 

dq° A F (ki + q) ■ ■ ■ A F (k r + q) = -in ^ — J[ S (k b + q 

a=l ZPa b+a 



-U)a±Pa 



-ITT 



-ITT 



duS^ - — [S(u + uj a +p a ) + 5(uj + uj a ~ Pa)) TT S (k b + q) 

/oo r 
duj^25((k a + q) 2 ) l[S (k b + q). 
-OO 1 i _/_ 



a=l b^a 

Using (|3.32p gives the result. In the case r = 1, the two summands in (|4.19p do not 
decay fast enough for closing the contours. However, one can compute the contour 
integrals using (|4.18j) . to again obtain the result. □ 

This lemma shows that in a closed loop, a product of r Feynman propagators may 
be replaced by the product of one factor Pq and r — 1 factors So, if we take the sum 
over all such combinations. The condition that the momenta k±, . . . , k r should all be 
different can be justified as follows: When computing scattering amplitudes, the mo- 
menta k\, . . . , k r are expressed in terms of the incoming and outgoing momenta. Then 
the momenta k a are are all different except at the poles of the scattering amplitudes. 
Such poles correspond to resonances; they need to be described non-perturbatively 
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with a resummation technique. If the loop under consideration is contained inside an- 
other loop, one also needs to integrate over one or several of the momenta k a . In this 
case, the poles of the Feynman diagrams if some of the momenta k±, . . . ,k r coincide 
are of no relevance, as only the divergences for large momenta are of interest. 

Thus, in order to decide whether the classical and the quantum loop diagrams agree 
quantitatively, it remains to study the combinatorics. Let us begin with the tadpole 
diagram (see Figure [1] (b)). It arises when integrating the simple tree diagram (see 
the left of Figure [5]) over S, if two of the incoming fields are the stochastic field £. As 
there are three possibilities for the choice of the two £-legs, and taking into account 
the factor 1/6 in (|3.29p . we get a total combinatorial factor of 1/2. Hence Lemma 14.31 
allows us to replace the factor Pq in the classical tadpole by 



As the factor (— i) is already taken care of in (I2.38p . we get complete agreement between 
the quantum and the classical tadpole if we choose (see (I4,5p and (I4.10p ) 



We next consider the simple loop diagram in Figure Q] (d). In the classical setting, 
this diagram arises by performing the stochastic integral of a second order tree diagram 
with two £-legs (see Figure [6] (b)). As the tree diagram has a combinatorial factor 
one (see Proposition 13. 5p . the inner lines are formed by the product X^SqPq. Due 
to the symmetry in the momenta of the two inner lines, Lemma 14.31 allows us to 
replace the inner lines by AfAf/(—2tti). Hence we again get complete agreement 
between the quantum and the classical diagram if we choose /3 according to (I4.20p . This 
argument applies just as well to simple loops inside larger diagrams, as is exemplified 
in Figure [6] (c). 

The situation changes with longer loops, as we now discuss in the example of the 
triangle loop in Figured] (d). The shown diagram is obtained in a unique way as the 
stochastic integral of a third order tree diagram with two £-legs. The combinatorial 
factor is again one. Similarly, one gets the two loop diagrams where the line Po is one 
of the other sides of the triangle. Thus we can again apply Lemma 14.31 In order to 
obtain complete agreement to the corresponding quantum diagram we need to choose 



We point out that this differs from (|4.20p by a factor of two. This shows that the 
expansion in terms of classical loop diagrams is not mathematically equivalent to 
quantum field theory. One difference is that longer loops have a weight which is 
double that of a tadpole or a simple loop. 

For higher loops (see Figure [6] (e) for the simplest example), one can apply induc- 
tively Lemma [4.3l (or a generalization where on the left factors Ap are replaced by So) 
to get a connection between the classical and the corresponding quantum diagrams. 
However, the combinatorial details of this analysis go beyond the scope of the present 
paper. 




1 




(4.21) 



4.4. The Zero-Point Energy of the Stochastic Field. A quantum field has an 
infinite zero-point energy, which is detectable via the Casimir effect. The analogous 
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quantity in our classical setting is the mean energy of the classical stochastic field £, 
given formally by 

'foO^E. (4.22) 



/ 



We point out that following our discussion in Section (|3.3p . this "zero point energy" 
is not directly accessible to an observer, as only energy differences can be measured. 
Nevertheless, it is conceivable that this energy could be measured in more realistic 
physical models which include fermions and thus make it possible to set up the Casimir 
experiment. With this in mind, it seems interesting to compute (|4,22p and to compare 
the result with the zero-point energy of a quantum field. 

The mean energy of the stochastic field (j4.22p is infinite for two reasons: First, the 
mean energy density is infinite in view of ()4.5p and (|4.10p , taking into account that the 
distribution Pq(x,u) is singular at x = y. Second, the integral of the energy density 
diverges because the spatial volume is infinite. In order to avoid the second problem, 
we consider the system as usual in a three-dimensional box of length L with periodic 
boundary conditions. Then the Fourier integral (|3.32p is to be replaced by the Fourier 
series 

fcS(27rZ/L) 3 1 1 

where we set co = \k\ and t = x° — y°. Using (I4.10P and (|3.2p . we obtain 

fce(27rZ/L)3 1 1 s,s'=±l 
fce(27rZ/L) 3 s=±1 fce(27rZ/L) 3 

Here the infinite sum clearly diverges. But to every momentum mode we can associate 
the mean energy 

E U 



e ik(x~y) 



27T/3 

Choosing f3 according to (|4,2ip gives perfect agreement with the zero point energy fko/2 
of a radiation mode of a quantum field as described by a quantum mechanical harmonic 
oscillator. 



5. Outlook 

The previous analysis shows that by considering the classical c/> 4 -theory in a sto- 
chastic background field, one gets a perturbative description which involves all the 
tree and loop diagrams of the corresponding quantum field theory. The main result of 
the paper is to demonstrate that loop diagrams appear in a purely classical context. 
By adapting our methods, this result immediately carries over to any other bosonic 
field theory in any space-time dimension. 

We also observed that for t/> 4 -theory, the combinatorics of the tadpole and of the 
simple loops differs from that of the longer loops (see Section 14. 3 fl . This result is quite 
sensitive on the theory under consideration. It might well be that for other theories, 
these combinatorial differences between the classical and the quantum loop diagrams 
disappear. In particular, it was crucial for our combinatorial considerations that there 
are distinct lines which begin at the same vertex i and both end at the same vertex j. 
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As such lines do not exist in QED, is conceivable that for QED there might be full 
agreement between classical and the quantum loop diagrams. 

In any case, working out the combinatorial details for a realistic model to higher 
order seems an important project for the future. Here we only make three remarks: 

(a) We note that for complex fields (like a complex Klein-Gordon field), the charge 
is an additional conserved quantity. This makes it possible to work in the clas- 
sical measurement process instead of the energy with the charge. An advantage 
of this procedure is that the charge is a Lorentz scalar, so that the time de- 
rivative in (|3.48p could be omitted. The reason why we preferred to work 
with a real field is that the components of the electromagnetic field are also 
real- valued. 

(b) We remark that in order to extend our methods, one could replace the Gauss- 
ian measure (|4.ip by a more general probability distribution on the classical 
solutions. This would change the detailed form of the higher order loop dia- 
grams. The non-Gaussian probability distribution would have to be justified 
as arising from the nonlinear classical dynamics prior to the scattering process. 

(c) We point out that all our methods only work for bosonic fields. The reason 
is that the Pauli exclusion principle and the transformation of half-integer 
spin fields under rotations seem to make it impossible to interpret the Dirac 
equation as a classical field equation. Thus in order to extend our methods to 
models involving fermions, one possibility is to couple a classical bosonic field 
to a second-quantized fermion field. An alternative is to describe the fermions 
in the framework of the fermionic projector [3J. 

We finally outline how our constructions could be helpful for the rigorous construc- 
tion of interacting quantum fields. To this end, let us assume that there is a physical 
theory for which there is complete agreement between the classical and the quan- 
tum Feynman diagrams. Then one could prove the existence of the quantum theory 
by rigorously constructing the classical theory in the presence of a stochastic back- 
ground field. More precisely, one could first introduce an ultraviolet regularization 
by replacing the covariance of the Gaussian measure (|4.10p by a regularized covari- 
ance P e / ' j3 (for example with a cutoff in momentum space; clearly, the regularization 
would violate Lorentz invariance). Then the initial values (j>o\t=-T (with cj)Q according 
to (|4.16p ) would be in a suitable Sobolev space. If one could prove global existence 
of solutions uniformly in £, integrating the solutions over H would give all interacting 
n-point functions. If one also succeeded in controlling the behavior of the classical 
solutions in the limit when the regularization is removed (this would probably require 
some kind of "non-perturbative renormalization method"), one would have proved 
non-perturbatively that the quantum field theory in the Lorentzian setting is well- 
defined and finite. In this way, it seems possible at least in principle to make the 
powerful methods of nonlinear partial differential equations applicable to one of the 
most outstanding problems of quantum field theory. 
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